Diffusive mass transfer by non equilibrium fluctuations: Fick's law revisited 
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Recent experimental and theoretical works have shown that giant fluctuations are present during 
diffusion in liquid systems. We use linearized fluctuating hydrodynamics to calculate the net mass 
transfer due to these non equilibrium fluctuations. Surprisingly the mass flow turns out to coincide 
with the usual Fick's one. The renormalization of the hydrodynamic equations allows us to quantify 
the gravitational modifications of the diffusion coefficient induced by the gravitational stabilization 
of long wavelength fluctuations. 
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Recent experiments have shown that giant non equilibrium concentration fluctuations are a universal feature of 
diffusion in binary liquid mixtures and macromolecular solutions || . The presence of long-range correlations in a 
diffusing liquid mixture is a surprising result, as diffusion is usually thought to give rise to a homogeneous mixing 
at the molecular level, other length scales playing no role in the process. The results of these experiments have 
been interpreted by means of an extension of Landau's Fluctuating Hydrodynamics || to time-dependent diffusion 
processes j|. 

Linearized hydrodynamics shows that the fluctuations are generated by thermal velocity fluctuations, due to the 
presence of a macroscopic non equilibrium concentration gradient. The mechanism is analogous to that first predicted 
|^|, H and then observed j?], ||, |j| [l(J for stationary thermal diffusion in a binary liquid mixture. Velocity fluctuations 
parallel to the macroscopic gradient displace parcels of fluid into regions with different concentration, thus giving rise 
to concentration fluctuations [ pd| . Since the fluctuations displace mass along the concentration gradient, one might 
wonder if they contribute significantly to the macroscopic diffusive mass flow. 

We will show that, quite unexpectedly, the net mass transfer due to non equilibrium fluctuations corresponds 
to the macroscopic Fick's mass flow. Traditionally the diffusive flow has been interpreted as due to the random 
hops of the molecules in the fluid. Here we show that Fick's flow is due to random hops of parcels of fluid whose 
dimensions span all the length scales from the molecular up to the macroscopic one. This is a surprising result, as 
usually thermal fluctuations are thought of as small perturbations of a macroscopic state, while here we show that 
the macroscopic flow is determined by a second-order effect in the fluctuations. The introduction of new degrees 
of freedom, the hydrodynamic fluctuations, determines an ultraviolet divergence of the diffusion coefficient, which 
we regularize by imposing a cut-off at large wave vectors. The regularization procedure involves the introduction 
of an arbitrary length scale. By imposing that the cut-off length scale is of the order of the correlation length of 
the mixture we find that the diffusion coefficient displays the usual Stokes-Einstein form and is of the same order 
of magnitude of the phenomenological one. However the introduction of an arbitrary length scale into the theory 
can be avoided by a renormalization of the hydrodynamic equations. The renormalization procedure allows us to 
quantify the modifications of the phenomenological diffusion coefficient due to the gravitational stabilization of long 
wavelength fluctuations. 

We will now derive the net contribution of the non equilibrium fluctuations to the mass transfer. The sample 
considered is a binary mixture of concentration c (x, t) under rest. A local fluctuation in the current of mass can be 
written as the second order product 

%(x,t) = p6c(x,t)6u(x,t), (f) 

where 5c(x, t) is the local concentration fluctuation and <5u(x, t) the local velocity fluctuation. 
The net contribution of the fluctuations to the diffusive current is obtained by averaging over the configurations of 
the system: 

j / (x,t) = p(5c(x ) i)5u(x,i)). (2) 

Equation (Q) shows that a current fluctuation takes place whenever a concentration fluctuation and a velocity 
fluctuation occur simultaneously at the same place. At thermodynamic equilibrium fluctuations in the velocity and 
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concentration are not correlated. Therefore the probability of a current fluctuation is small and with random sign, 
and the net flow j f averages to zero. Under non equilibrium conditions, however, velocity fluctuations are correlated 
to concentration fluctuations [Q 0, |(| . The presence of a macroscopic concentration gradient determines the existence 
of velocity- induced concentration fluctuations: a velocity fluctuation determines a small vortex inside the fluid and 
this brings parcels of fluid with a certain concentration into layers of fluid with a different concentration. On the 
contrary, the gravity force determines the opposite mechanism: whenever a concentration fluctuation is produced, the 
buoyancy force tends to displace the fluctuation. Therefore terms in the average of Eq. (||) are correlated and this 
determines a non vanishing mass flow j f . 

We will shortly derive an explicit expression for j/ by determining the cross-correlation properties of concentration 
and velocity fluctuations in the presence of a macroscopic concentration gradient. This task can be accomplished by 
using Landau's linearized hydrodynamics H. Basically one writes down the relevant hydrodynamic equations for a 
binary mixture. The equations are then linearized for small fluctuations of the hydrodynamic variables around their 
macroscopic values and random source terms are added to them. 

The equations thus obtained for a binary mixture where a uniform concentration gradient is present are B 

\ 5c (q, lj) (iuj + D q 2 ) = -5u (q, ui) ■ Vc + iq • F (q, ui) 

(3) 



\ <*u (q, w) (iw + vq 2 ) = /35c (q, w) (g - q*fj + iq A S (q, w) ' 

where Dq is the bare diffusion coefficient, v the kinematic viscosity, g the gravitational acceleration, (3 = p^ 1 (dp/dc) 
and c (q, lo) is the space-time Fourier transform for the concentration fluctuations defined by 

5c (q, oj) = - J-j f ■ x + lujt 5c (x, t) dxdt, (4) 
(27r) J 

and an analogous expression holds for the velocity fluctuations. 

In Eq. (^j)F(q, w) and S(q, u>) represent random source currents for concentration fluctuations and velocity fluctu- 
ations, respectively. As customary, their correlation properties are assumed to retain their equilibrium values @: 

(F i (q,o;)F J (q',u;'))=<Si J <S(q + q')^ + ^)^ffy (j^j ^ (5) 
(Sitq^-fa'V)) ^^^q + qO^ + c')^^, (6) 



(F i {q,oj)S j (q!,u;'))=0. (7) 

Let us briefly comment Eq. (|^). The left terms specify that concentration and velocity fluctuations decay with 
lifetimes 1/Doq 2 and 1/vq 2 , respectively. The right terms are source terms for the fluctuations. Beyond the random 
currents F(q, lj) and S(q,oj), they contain the cross-coupling between velocity and concentration fluctuations. The 
cross-coupling term in the concentration equations shows that a velocity fluctuation 5u couples with the macroscopic 
concentration gradient Vc, thus giving rise to a concentration fluctuation. The cross coupling term in the velocity 
equation represents the acceleration acting on a concentration fluctuation due to the buoyancy force. We have 
assumed that the fluid is incompressible; therefore the acceleration acting on the fluctuation is the component of 
gravity perpendicular to the wave vector q. 

The explicitation of 5c and 5u in Eq. (H) yields 

!(iw + ^ 2 )q-F(q,w)-hAS(q,w)-Vc 
5c (q, cj) = r (8) 



(ilo + D Q q 2 ) [iu> + vq 2 ) + 13 (g - qaf ) ■ Vc 



and 
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Su(q,u) = 



*q-F(q,w)0(g-qaSP) 


+ i(w + fl 9 2 )qAS(q, W ) 


(iw + £) g 2 )(^ + ^ 2 )+/3 


(g - q*#) 


• Vc 
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We are now in the position to evaluate the mass flow j/. In the reciprocal space Eq. (||) becomes 

if = p J (Sc(q, u)Su(q', a/)) dqdq'dwda/. 



(9) 



(10) 



The integrand in Eq. ( |10| ) can be decomposed as the sum of terms containing correlations of the random currents 
F and S. Due to (0) the only non vanishing terms are those originated from self-correlations: 



(<5c( q ,^u(q>')> = (<5c(q,^u(q>')}f + (Mq^)Mq>')> 



S ' 



where 
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(6c(q, w)iu(q', u')) F = S(q+q')6(u + u') *£^5° ( ?f) (3 



q 2 (iu + vq 2 ) (g - q^3 



(iu + D q 2 ) (iu + vq 2 ) +fi(g- q^F J ■ Vc 



and 
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(<5c(q, w)5u(q', w'))s = (q + q') 5 ( w + w ') 



<z 2 (^ + ^ 2 ) (Vc-q^ 



Vc q 
T~ 



8tt 4 p 



(iw + D a q 2 ) (iu + vq 2 ) + (3 (g - q^j ■ V( 



(13) 



The denominator of Eqs. jl^-p^|) is an even function of u; therefore in the evaluation of the integral in Eq. ( |l0| ) 
we can neglect odd functions of u in the numerator. 

By combining Eqs. (]l0|— p~3|) and by assuming that z = g/g = Vc/|Vc|, the current jj due to non equilibrium 
fluctuations is finally determined: 



where 



J/ = -P D f 



(14) 



D, 



K b TvDq 
8tt 4 p 



^(--qf 1 



(iu + D q 2 ) (iu + vq 2 ) + figVc I i — q 



r dqdaj. 



(15) 



Surprisingly the flow j/ displays the customary dependence of Fick's flow on the driving force Vc and on the 
barodiffusion flow. 



Equation ( |15| ) can be rewritten by factorizing its denominator as (u + uj_) (u + u_) and by approximating the 

two roots as u\ — D^q A |l + [q/qao (<?)] 4 (z — q^J ' ^} ano - w - = i/2 <? 4 - Basically this approximation holds true 
for v >> D, a condition fulfilled by most binary liquid mixtures; further details are discussed in Ref. B: 
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Df 



where 



K B T 



pv 



q 4 + q RO {gf 



-dq, 



(16) 



quo (g) = 



g(3Vc 



vD n 



(17) 



is a roll-off wave vector due to the gravity force || . As thoroughly discussed in Ref. Q , the roll-off wave vector 
marks the transition from the diffusive relaxation of fluctuations at large wave vectors to the gravitational stabilization 
of fluctuations at smaller ones. At the roll-off wave vector the diffusional and gravitational fluctuation relaxation times 
coincide. 



By introducing polar coordinates Eq. (16) becomes 



(2ir) pv Jo Jo Jo < 



q 4 sin 3 (0) 



qno(g) sin (0) 



(18) 



The integral on q is divergent. Therefore we have regularized it by introducing a brute- force cut-off at the arbitrary 
wave vector Q. The evaluation of the integrals in tp and q yields 
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By assuming that q R o (g) /Q *C 1 Eq. (|19|) becomes 



\f2-Kq RO r^—r ( q RO Vsm9 
1 : — Vsint' + O 



(19) 



D ? (9) = - ^TT^-mo (g) + QO (q RO (g) IQ) 4 , 



3 {2^Y pv 



where 



4 (2tt) pv 



(20) 



I = 



d9sin 7/2 9 re 1.248599 



In the absence of gravity q R o = 0, and Eq.(p20|) becomes 



D, 



2 K B T 
3 TrpvA ' 



(21) 



(22) 



where A = 2ir/Q can be roughly assumed to be of the order of correlation length of the fluid. It should be noticed 
that the diffusion coefficient Df of Eq. (^2|) does not depend on the bare mass diffusion coefficient Do which appeared 
in the linearized fluctuating hydrodynamic equations (||) and (||). 

Quite remarkably, Eq. ( ]2^ ) exhibits the usual Stokes-Einstein dependence on temperature, viscosity and correlation 
length, the proportionality factor being of the right order of magnitude. This result, combined with Eq.(|l4|), strongly 
suggests that the mass transfer due to fluctuations and Fick's mass flow can be identified. This is a striking result, 
as it shows for the first time that the diffusive mass transfer is entirely determined by non equilibrium fluctuations 
occurring at length scales ranging from the microscopic up to the macroscopic one. 

Although it is physically reasonable to assume that wave vectors larger than a typical molecular one have no 
meaning, the imposition of a cut-off involves the introduction of an arbitrary length scale, which in turn requires to 
introduce a microscopic model for the molecular interactions. This necessity can be eliminated by a renormalization of 
the linearized hydrodynamic equations. The normalization procedure we will discuss in the following closely mirrors 
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the renormalization of mass of Quantum Field Theory and of Classical Electrodynamics |Tj] . The overall diffusive 
mass flow is the superposition of the contribution due to fluctuations j / and the bare fickean mass flow: 



J (.9) = (9) + J? (<?) = -P (Df (g) + D%) 



V »'g ( |^ 

<9/i/ p T 



The renormalization condition is 



(23) 



j (.<? = 0) = -pDVc, 



(24) 



where D is the phcnomcnological diffusion coefficient in the absence of gravity. 

By combining Eq. (E3h and Eq. (J24|) the renormalization condition can be rewritten as 



D$ = D- D Q S ( 5 = 0), 
and, by inserting Eqs. (EG,E5h into Eq. (E3l), the renormalized mass current is 



(25) 



5r (d) = ~pD R (g) 



Vc-/3g 



dc 



(26) 



where Dr is the renormalized diffusion coefficient 



D R (g) =D- V2jl B q RO (g) . (27) 

In Eq. ( ]27| ) the Q dependence has been eliminated by letting Q — > 00, because Dr does not diverge. The microscopic 
details have been reabsorbed into D, the only experimental parameter needed by the theory in order to obtain the 
value of the diffusion coefficient in the presence of gravity. 

Equation (5?]) shows that the diffusion coefficient is smaller due to the presence of gravity, the additional term 
having the Stokes-Einstein form, the only length scale involved being 2ir jqRo- 

This modification of the diffusion coefficient is due to the fact that, in the presence of gravity, fluctuations with 
wave vector smaller than qro are stabilized and therefore do not contribute to the mass transfer. Therefore qRo 
acts as an infrared cut-off for the fluctuations. As the length scales stabilized by gravity are large, the relative 
correction introduced by gravity is very small, being of the order of qRo£,, where we have assumed a Stokes-Einstein 
dependence of the phenomenological diffusion coefficient D on the correlation length £ of the fluid. For an ordinary 
liquid mixture qRo can be as high as 10 3 cm -1 and £ is of the order of 10~ 7 cm. Therefore the gravitational correction 
is estimated to be roughly one part in 10 4 . In the case of critical fluids and macromolecular solutions the correlation 
length can grow very large, thus suggesting that the gravitational correction to the diffusion coefficient can get large. 
However, in order for our hydrodynamic description to be valid, the correlation length has to be significantly smaller 
than 1 j qRo- Therefore the largest correction to the diffusion coefficient reliably predicted by this description can be 
roughly estimated to be of the order of 10%. Such a variation in the diffusion coefficient should be observable by 
performing microgravity experiments on critical fluids. 
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